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We study the trajectories followed by a particle subjected to weak noise when escaping from the 
domain of attraction of a stable fixed point. If detailed balance is absent, a focus may occur along 
the most probable exit path, leading to a breakdown of symmetry (if present). The exit trajectory 
bifurcates, and the exit location distribution may become 'skewed' (non-Gaussian). The weak-noise 
asymptotics of the mean escape time are strongly affected. Our methods extend to the study of 
skewed exit location distributions in stochastic models without symmetry. 
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A particle moving in a force field, but weakly perturbed 
by external noise, will spend most of its time near sta- 
tic', ble fixed points of the force field. But the particle will 
^ I occasionally undergo a large fluctuation: it will leave the 
' . basin of attraction of one such point, and enter that of 
00 ' another. The timescale on which such fluctuations occur 

■ grows exponentially in the low-noise limit. 

I I In this limit, the particle typically follows a unique tra- 
^ . jectory in its final, successful escape attempt; transverse 
' deviations from this path become increasingly unlikely 
• as the noise strength goes to zero. This trajectory is 
I referred to as the most probable escape path (MPEP), 
. and its properties govern the asymptotic behavior of the 
mean escape time. In our earlier work we suggested 
0^ the possibility of a focusing singularity: as one moves 
out along the MPEP, transverse fluctuations become in- 
creasingly less stable, leading to a breakdown of stability 
before the boundary of the basin is reached. This hith- 
Q . erto unexplored phenomenon has profound consequences 
for the particle's escape behaviour. It frequently occurs 
when the force field is nonconservative, or in general 
J-^ ■ when detailed balance is absent. Its wide prevalence has 
^ . not previously been recognized. 

In this Letter we study a simple example to demon- 
. strate the effects of focusing singularities on exit phe- 

■ nomena. We show that when such a singularity occurs, 
the MPEP bifurcates (in the simplest case), leading to 
a type of broken symmetry. This transition affects not 
only the prefactor of the mean first passage time (MFPT) 
to the boundary, but also, its exponential dependence 
on the noise. It also induces an unusual 'skewed' (non- 
Gaussian) distribution of exit points on the boundary. 

For simplicity we restrict ourselves to the case of an 
overdamped classical point particle moving in two dimen- 
sions and subject to additive isotropic white noise w(t). 
In the overdamped case the deterministic forces acting 
on the particle are described by a drift field u(x), so that 

dx,{t) ^Ui{-x{t))dt + e^/^dw^{t), i = l,2 (1) 

with e the noise strength. The corresponding Fokker- 
Planck equation for the probability density is 



1 



p=(e/2)VV-V-(pu). 



(2) 



We shall employ the technique of matched asymptotic ex- 
pansions described in Refs. 0J^. We impose absorbing 
boundary conditions on the 'separatrix' (the boundary 
of the basin) and study the behavior, as e — > 0, of the 
quasistationary density pi{x,y), i.e., the slowest decay- 
ing eigenmodc of the Fokkcr-Planck operator. In the low- 
noise limit, the properties of pi can be used to calculate 
both the MFPT and the distribution of exit points ^ . 

We specialize to the case of a drift field with a stable 
fixed point at (a;5,0), whose basin of attraction is the 
entire right-half plane. We assume symmetry about the 
X-axis, and that there is only one fixed point on the y- 
axis: a saddle point at (0,0). A simple drift field with 
this structure is 

\i{x,y) ^ {x ~ ~ axy^ , - py ~ x^y). (3) 

This force field is conservative, i.e., is the gradient of 
a potential, only if a = 1. Nonetheless for any p > 
the field conforms to the assumptions, with 2:5 = 1. The 
symmetry of u suggests that the MPEP, which must con- 
nect {xs, 0) and the saddle, lies along the x-axis, and for 
small a > that is the case The case when exit 
occurs over an unstable fixed point is also of interest; 
additional phenomena may emerge there because focus- 
ing on the separatrix (rather than the MPEP) can occur. 
We defer that case to a later paper. 

Away from both the stable point and the separatrix, 
pi can be approximated by a WKB form [^|j3) 

pi{x, y) - K{x, y) exp {~W{x, y)/e) . (4) 

K[x, y) satisfies a transport equation, and W{x, y) an 
eikonal (Hamilton- Jacobi) equation: i/(x, VW) — 0, 
with -ff (x, p) = ip^ -I- u(x) • p the Wentzell-Freidlin 
Hamiltonian [Q. So for any point {x,y) in the basin of 
attraction of {xs, 0), W{x, y) will be the action of a zero- 
energy classical trajectory, governed by this Hamiltonian, 
extending from {xs, 0) to (x, y). In general, computation 
of W{x,y) will require a minimization over the set of 
zero-energy trajectories terminating at (x,y). MPEPs 
are accordingly the zero-energy trajectories from (xs,Q) 
to (0, 0) of least action. The MPEP action is the expo- 
nential growth rate, in the low-noise limit, of the MFPT. 

For non-gradient drift fields u, it will be a common 
occurrence for two or more classical trajectories ema- 
nating from (a;s,0) to cross at some point. The locus 
of these points (equivalently, the set of points reachable 
from {xs, 0) via more than one zero-energy trajectory) is 
known as a caustic |^-|^ . A focus is a critical point from 
which a caustic emanates. Pictorial examples of this are 
given in Refs. j^J^,^. In these examples the caustics are 
found far from fixed points, MPEPs, or separatrices and 
have little effect on exit phenomena. However they may 
constitute regions of anomalous probability density, since 
the WKB expansion breaks down there . 
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Our concern here is what happens when a focus ap- 
pears on the MPEP itself, making the action W non- 
unique at foUowing points. Because the MPEP is initially 
assumed to lie along the x-axis, we study the validity of 
the WKB approximation by introducing the expansion 

W{x,y)^fo{x)+f2{x)y^ + oiy^). (5) 

Physically, f2{x)-'^^^e^^^ is the transverse lengthscale 
within which pi{x,y) is nonnegligible. f2{x) satisfies a 
nonlinear Riccati equation, described in Refs. UM- 

We begin by considering the drift field of (3^ when 
jjL — 1. For all a < 4, we find numerically that remains 
positive and finite from a; = 1 to a; = 0, converging to 1 
as the origin is approached. (Because the origin is a fixed 
point of u, it is formally reached in infinite time; see the 
discussion in Rcf. [^). However when a = 4, /2 con- 
verges to zero, signalling the appearance of a transverse 
'soft mode,' or instability. (This is accompanied by a 
divergence of the prefactor K in (^).) For all a > 4, 
/2 reaches zero at positive x, and then plunges to —oo in 
finite time (hence at strictly positive x) . This signals the 
presence of a singularity in W . 

Fig. |l|(a) displays this graphically. As one moves 
from (1,0) toward the origin, the iso-action curves be- 
come non-convex at the point where /2 goes negative. 
They become increasingly pinched near the x-axis, lead- 
ing to a cusp singularity at (a;_F,0), where /2 diverges. 
Beyond this point W remains continuous, but is not con- 
tinuously differentiable across the x-axis. It must be com- 
puted by minimizing over trajectories, and different tra- 
jectories are employed as y — > 0+ and ?/ — > 0~. 

Fig. |l|(b) makes this clear; it shows the zero-energy 
classical trajectories, derived from the Wentzell-Preidlin 
Hamiltonian by integrating Hamilton's equations. As x 
is decreased from 1, off-axis trajectories eventually re- 
converge on the axis, signalling that {xf,Q) is a focus. 
(It is not hard to show this analytically, but we omit 
the proof here.) Every point (x,0) with x < xp, includ- 
ing the saddle point at the origin, is reachable by three 
zero-energy trajectories: the original x-axis MPEP, and 
two symmetrically placed off-axis trajectories, with lesser 
(degenerate) actions. 

The dependence of the focus xp ona, for the particular 
drift field (|^), is of some interest. When a = A, xp = 0. 
As a increases beyond 4, x^? increases continuously from 
zero but never reaches xs- When a = 9, a new focus 
appears at the origin so that two foci occur on the x-axis. 
Both of these move to larger a; as a increases. We have 
found numerically that a new focus appears at a — n^, 
n = 2,3,4, .. ., but have not yet explored in detail the 
set of zero-energy trajectories when multiple foci occur. 
We restrict ourselves below to the case of a single focus, 
i.e., if /i = 1 we consider only 4<a<9. If/iy^l the 
range of allowed values of a will differ. 

Upon encountering a single focus, the MPEP bifur- 
cates. The path of minimum classical action — hence. 
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the MPEP — to any point {x, 0) with x > xp from (xg, 0) 
remains the x-axis. For x < xp, however, the former 
MPEP (i.e., the x-axis) remains an extremum of the ac- 
tion but is no longer a local minimum Q . There are two 
new MPEPs, related hy y ^ —y. The symmetry has 
been broken: the drift field and the equations of motion 
are symmetric about the x-axis, but each MPEP is not. 

The new (curved) MPEPs, and their common action, 
are easily computed numerically. Since their action is 
the exponential growth rate of the MFPT, the bifurca- 
tion greatly affects the low-noise asymptotic behavior. 
In Fig. 1^ we plot the classical action of the MPEP(s) be- 
tween the stable and hyperbolic points (and therefore the 
asymptotic slope of the MFPT vs. 1/e on a log- log plot) 
vs. a for the case ^ — 1 in the drift field (§). The onset 
of focusing at a = 4 is evident. The MFPT prefactor, 
one can show, is also singular ai a — A. 

To study the effects of focusing on the exit location 
distribution near the saddle, we set up a covariant for- 
malism for computing the WKB behavior of pi along a 
curved MPEP. The Hamilton- Jacobi equation is 

= Hix\p,) = ^D^'p^Pj + u\^)p„ (6) 

where pi = dW/ dx^ , an index is summed over when re- 
peated, and the diffusion tensor D may be anisotropic 
(but we consider only isotropic D here). We adopt the 
notation pij = dpi/dx^; note that W,ij = pij. By twice 
differentiating the Hamilton-Jacobi equation we arrive at 
the evolution equation for the second derivatives of W 
along any zero-energy trajectory: 

= -D'''W,mW,i, - u\, W^kj - u\j W,k^ - u'.,^pi 

(7) 

where the dot signifies a derivative with respect to transit 
time along the trajectory. 

Eq. (|^) describes the evolution of fluctuations in all 
directions of the action along the MPEP, and represents a 
generalization of our previously derived nonlinear Riccati 
equation (to which it reduces when i — j — y and 
the MPEP is along the x-axis). Similarly, the covariant 
generalization of our transport equation for A" is 

k/K ^-u\^ \D'i W^^j . (8) 

It is amusing to note that the diffusion tensor plays a role 
similar to that of a metric tensor, with additive noise 
(constant D) corresponding to flat space, while multi- 
plicative noise introduces an effective curvature 

In the vicinity of the stable and hyperbolic fixed points 
{t — > — oo,oo above) the momentum and therefore the 
inhomogeneous term of (^ vanish. We see that in these 
limits we must solve an algebraic Riccati equation 

D'^'W^mW^ij + u\, W^k, + u^J W,k^ - 0. (9) 
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Eq. was also derived by Ludwig and Mangel in the 
context of stochastic ecology The quasistationary 

density in the far field of the 'diffusive' lengthscale, of 
size 0(e^^^), surrounding each fixed point is accordingly 
(up to a constant factor) 

pi(x) ~ cxp [-(x - xo)C-i(x - xo)/2e] (10) 

where xq = (x5,0) or (0,0), and the inverse covariance 
matrix C^^ — {C~^)ij — W^ij is determined by the ma- 
trix equation (^. C^^ will be positive definite at the 
stable point, but not at the saddle point. 

Our procedure is to compute W(x) and K{x.) by (nu- 
merically) integrating the coupled equations (0) and (H) 
along the curved MPEP. Close to the two fixed points, 
we match to a quasistationary density which (a) is a so- 
lution of the Fokkcr-Planck equation, and (b) matches 
with (^ on the ©(e^/^) lengthscale. 

We find that the asymptotic behavior of exit phenom- 
ena falls into three qualitatively different categories de- 
pending on whether the parameter /i in the drift field is 
greater than, less than, or equal to 1. This is because 
in the vicinity of the saddle at the origin, u « (x, —fJ.y). 
In this linearized field Hamilton's equations are easily 
solved. (All trajectories possess a reflection symmetry 
about the x-axis; when we refer to a trajectory in the first 
quadrant, note that it has a mirror image in the fourth.) 
The zero-energy trajectory which approaches the origin 
from above follows a curve y cx x''; hence the MPEP 
asymptotically approaches the x-axis for fi > 1, the y- 
axis for < 1, and comes in at a fixed angle for /i = 1. 
[Cf. Fig. |I](b).] We consider these cases in turn. 

> 1: Very near the saddle point at (0,0), this 
case resembles the case of unbroken symmetry described 
in Ref. ||^. The action there can be approximated as 
W{x,y) « —x^ + ^y^. The prefactor K simply tends to 
a nonzero constant. Because the MPEPs approach the 
origin asymptotically close to the x-axis, matching in the 
vicinity of the saddle remains the same as in the unbro- 
ken case. Up to a constant factor, the quasistationary 
density pi on the x,y = 0(e^/^) lengthscale wiU be 

exp(xV2e)y2(i2:v/2A)] exp(-/V/e) (H) 

where 2/2(^1 •) is the odd parabolic cylinder function 



of index i. This agrees with (10) in the far field, and 
falls to zero as x ^ due to absorption of probabil- 
ity. The exit location distribution will be asymptotically 
Gaussian, with standard deviation e^/^/-y2/i, as in the 
unbroken case. We stress, however, that in this and the 
following two cases the exponential dependence of the 
MFPT on the noise depends on proper identification of 
the a-dependent off- axis MPEPs, as does the prefactor. 

/i < 1: Here the MPEP approaches the origin asymp- 
totically along the y-axis. However, the momentum of 
the MPEP asymptotically aligns itself along the x-axis. 
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(This is true in all three cases.) This can easily be seen by 
solving Hamilton's equations of motion near the origin. 
Because the quadratic fluctuations of the action around 
the MPEP are perpendicular to the momentum, we have 
the unusual situation where 'transverse' fluctuations of 
the escaping particle around the MPEP become asymp- 
totically aligned with the MPEP itself: they become lon- 
gitudinal. The quasistationary density and exit location 
distribution are, moreover, spread over an anomalously 
large region. Absorption takes place in a boundary layer 
of width 0(e^/^), but because the MPEP approaches the 
origin along a curve y (x the appropriate lengthscale 
in the y-direction is y = 0(e^/^), not y = 0(e^/^). 

This becomes clear when W{x,y) and K{x,y) are ap- 
proximated near the origin. Integrating the evolution 
equation (^, in the linear approximation, yields 

W{x, y) « -2x{\y\/Af'^ + M/Af'^ (12) 

where A is the limit of the ratio y/x^^ along the MPEP 
as it approaches the origin. Since /i < 1, this action is 
of higher order than quadratic in the displacements x, y. 
Equivalently, the inverse covariance matrix C^^ of ( |l0| ) 
vanishes. This vanishing is always allowed \WAj = 
is a solution of the algebraic Riccati equation (^)] but 
it occurs only when /.t < 1, leading to an anomalously 
large lengthscale. The behavior of the prefactor K may 
be computed from (^); along the MPEP, it is asymp- 
totically proportional to x^"'^. This is an unusual sit- 
uation: strictly speaking, the 'frequency factor' Pp2| 
K{0, 0) equals zero, yet all physical exit quantities re- 
main nonzero and finite. The MFPT prefactor, however, 
turns out not to be anomalous; it is independent of e. 

To compute the exit location density and the MFPT 
asymptotics, an expression for the quasistationary den- 
sity pi on the X — 0(e^/^), y — 0(e^/^) lengthscale is 
needed. It is easily seen to be (up to a constant factor) 

sinh[2x(|2;|M)i/M/e] exp[-{\y\/Af/'^/e] 

which both solves the Fokker-Planck equation to leading 
order, and matches up with the WKB approximation (^) 
when W is given by (|l|) and K ~ x'^''' along the MPEP. 
Since the exit location probability density is proportional 
to dnPi, i.e., dpi/dx{x ~ 0), this gives an asymptotic 
exit location distribution f{y) dy on the y-axis, with 

/(2/)(x|y|(2/^')-iexp[-(|y|M)2/Ve] (13) 

So on the 0(e^/^) lengthscale to either side of the saddle, 
the exit location has a bimodal (symmetrized) Weibull 
distribution ||l^, with shape parameter 2/ p. (Equiv- 
alently, Ij/P^'' at exit time becomes exponentially dis- 
tributed in the low-noise limit.) This 'skewing' phe- 
nomenon (the particle tending to exit to one side of a 
saddle, in a non-Gaussian way, rather than at the saddle 



itself) was discovered by Bobrovsky and Schuss |14|, and 
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has been further investigated ll^Jl^]. However, in our 
simple model we are able to work out the skewed exit 
location distribution exactly. Of course the scaling pa- 
rameter A = A(a, fi) must be determined numerically, 
by integrating Hamilton's equations. The physical inter- 
pretation of the exponential distribution remains unclear. 

^ — 1: As in the /i > 1 case, absorption occurs on 
the x,y = 0{e^^'^) lengthscale. But a new feature enters: 
the inverse covariance matrix C^^ is no longer uniquely 
determined by (^). There is a free parameter, which 
must be computed numerically by matching. The asymp- 
totic expression for the exit location distribution involves 
parabolic cylinder functions. We defer the details. 

We have dealt only with a drift field symmetric about 
some axis. Our results (excepting those relating to bro- 
ken symmetry) and methods are more general and apply 
to drift fields without such symmetry. So they should 
apply, in particular, to the asymmetric models (with 
anisotropic D) reviewed by Bobrovsky and Zeitouni 
For such models the MPEP will typically be curved even 
in the absence of focusing, and skewed exit distributions 
will frequently occur. But the presence of a focus will 
still be signalled by a breakdown in transverse fluctu- 
ations, at which point qualitatively new behavior will 
emerge. One amusing prospect is to study drift fields 
with a 'symmetry parameter.' As this is tuned through 
zero one might observe a jump in the exponential growth 
rate of the MFPT, much as a jump in magnetization re- 
sults when an applied field passes through zero. 
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FIG. 1. (a) The curves of constant action W surrounding 
the stable fixed point {xs,0) = (1,0), revealing the presence 
of a singularity in W at {xf,0) ~ (0.3,0). Here n — 1 and 
a — 5. (b) The classical trajectories of zero energy, emanating 
from (1,0), from which the iso-action curves are computed. 
The two MPEPs are the off-axis trajectories incident on (0, 0). 

FIG. 2. The action of the MPEP(s) for the drift field (||), 
as a function of the parameter a when fi = 1. The MPEP 
extends along the x-axis only if a < 4, and no focus is present. 
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